Correlated two-photon scattering in cavity optomechanics 
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We present an exact analytical solution of the two-photon scattering in a cavity optomechani- 
cal system. This is achieved by solving the quantum dynamics of the total system, including the 
optomechanical cavity and the cavity-field environment, with the Laplace transform method. The 
long-time solution reveals detailed physical processes involved as well as the corresponding reso- 
nant photon frequencies. We characterize the photon correlation induced in the scattering process 
by calculating the two-photon joint spectrum of the long-time state. Clear evidence for photon 
frequency anti-correlation can be observed in the joint spectrum. In addition, we calculate the 
equal-time second-order correlation function of the cavity photons. The results show that the ra- 
diation pressure coupling can induce photon blockade effect, which is strongly modulated by the 
phonon sideband resonance. In particular, we obtain an explicit expression of optomechanical cou- 
pling strength determining these sideband modulation peaks based on the two-photon resonance 
condition. 

PACS numbers: 42.50.Pq, 42.50.Wk, 42.50.Ar, 07.10.Cm 
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I. INTRODUCTION 

The realization of photon-photon interaction at few- 
photon level has been a research subject of major interest 
in quantum optics 0, . The significance of few-photon 
interaction exists not only for studying the foundations 
of quantum theory, but also for applications in quantum 
information science. Specifically, an important goal in 
the current field of research is the control of two-photon 
correlations. Such a problem has been discussed in var- 
ious nonlinear systems [3H9j. For example, it has been 
shown that photon-photon interaction can be achieved 
in nonlinear Kerr media, and interesting quantum corre- 
lation effects such as photon blockade [l|, 0, LL(| appear in 
the strong-coupling regime. In particular, a strong Kerr 
nonlinearity can be achieved through the interaction be- 
tween light and atoms [l|, HJ . 

Interestingly, an optomechanical cavity [ll| - [l3| driven 
by radiation-pressure can be mapped to a problem with 
a Kerr- type interaction (l34l6| , and hence provides a dif- 
ferent class of systems to explore phenomena of interact- 
ing photons through the control of mechanical motion 
of the cavity mirrors. Indeed, Rabl pjj has examined 
the photon blockade effect in a continuously driven cav- 
ity optomechanical system operated in the single-photon 
strong-coupling regime. In such a regime, a single cavity 
photon can significantly change the resonant frequency of 
the cavity field and leads to a range of interesting effects, 
such as quantum state preparation [3, [l5|, [Til. Iigl , mul- 
tiple mechanical sidebands [2(|, scattering |2l|, single- 
photon cooling [221 ] , and optomechanical instability [23j j . 
We note that several recent experimental systems in cav- 
ity optomechanics [2^ - |26| are approaching the single- 
photon strong-coupling regime. 

In this paper, we investigate the quantum dynamics 
of a photon pair interacting with a moving mirror in a 



cavity. Different from previous studies (171 . |20| in which a 
continuously driven system is treated perturbatively, our 
system confined in the two-photon subspace is exactly 
solvable and therefore provides a fundamental configu- 
ration to study how two photons can become correlated 
via optomechanical coupling. We shall focus on the scat- 
tering problem in which a two-photon wave packet is in- 
jected into the cavity. By treating the cavity field, the 
moving mirror, and the field continuum outside the cav- 
ity as a whole system, we can solve the quantum state 
evolution via the Wigner-Weisskopf method. From the 
long-time solution, we identify quantitatively the proba- 
bility amplitudes associated with various scattering pro- 
cesses inside the cavity. 

We emphasize that our exact solution provides a com- 
plete description of the full quantum state, including the 
system's environment, and this is not directly attainable 
in previous related studies based on the master equation 
of the system's reduced density matrix f2l\ or approx- 
imated operator solutions in Heisenberg's picture [l7| . 
The knowledge of the full quantum state enables us to 
learn about the details of underlying physical processes 
as well as the quantum structure of interacting photons 
in optomechanical cavities. In this paper, we also cal- 
culate the joint spectrum of the two scattered photons. 
From the joint spectrum, we can see clear evidence for the 
photon frequency anti-correlation, which is a qualitative 
evaluation of the two-photon correlation. In addition, we 
calculate the equal-time second-order correlation func- 
tion of the cavity photons at transient times as a quanti- 
tative measurement of the two-photon correlation. The 
results indicate that photon blockade is strongly modu- 
lated by the optomechanical coupling strength as well as 
the Franck-Condon factors involved. More importantly, 
we find out the explicit modulation rule of the optome- 
chanical coupling strength based on the two-photon res- 
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FIG. 1: (Color online) (a) Schematic diagram of a Fabry- 
Perot-type optomechanical cavity formed by a fixed end mir- 
ror and a moving end mirror, (b) The energy-level structure 
(unsealed) of the optomechanical cavity (limited in the zero-, 
one-, and two-photon subspaces). 
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II. THE MODEL 

The system under consideration is a Fabry-Perot-type 
optomechanical cavity formed by a fixed end mirror and 
a moving end mirror [see Fig. Ufa)]. We focus on a 
single-mode cavity field, which is coupled to the mechan- 
ical oscillation of the moving mirror via radiation pres- 
sure. The moving end mirror is assumed to be perfect 
and the fixed one is partially transparent. In a rotating 
frame defined by the unitary transformation e ~ 4jv " c * with 
N = a^a + Jq° c^Ckdk being the total photon number op- 
erator, the Hamiltonian of the whole system including 
the optomechanical cavity and the environment reads as 

/>oo 

Hj = fruiM^b — hgoa^a(b^ + b) + / HAkC k Ckdk 

Jo 

/>OC 

+h£ / (c|.a + a ] c k )dh. (1) 
Jo 

Here a (cr) and b (b') are, respectively, the annihilation 
(creation) operators of the cavity field and the moving 
mirror, with the respective frequencies uj c and ujm- Ck 
(c k ) is the annihilation (creation) operator of the continu- 
ous field mode k outside the cavity with the resonant fre- 



quency Wfe (Afe = uj k —io c is the detuning). The radiation- 
pressure coupling appearing in the second term of Hi is 
characterized by the coupling strength g , and the cou- 
pling between the cavity field and the outside fields is 
modeled by the hopping interaction with the coupling 
strength £ = ^J^ c /2t: is the cavity- field decay rate). 
Since the decay rate of the mechanical resonator 7m can 
be much smaller than j c , we will neglect the dissipation 
of the mechanical resonator in our discussions. This is 
justified as long as the scattering processes are completed 
in a time much shorter than 7m • 

In this paper, we will focus on the two-photon scatter- 
ing problem. Because the total photon number is a con- 
served quantity, the Hilbert space for the photon part in 
this problem is spanned by three types of basis vectors: 
|2)„|0), |l>a|l*>, and |0>„|l p , 1,), where |2) a |0) stands for 
the state with two photons in the cavity and the outside 
fields are in a vacuum, |l) |lfc) represents the state with 
one photon in the cavity, one photon in the k mode of 
the outside fields, and the state |0) a |l p ,l g ) means that 
there is no photon in the cavity, two photons are in the 
p and q modes of the outside fields. By using these basis 
vectors for the fields, the state vector of the total system 
at time t is denoted by 

|$(<)) = f> m (i)|2) Q |0)|m(2)} & 

+ V / dfcB mifc (t)|l)„|l fc )|ro(l)) 6 
m=0 Jo 

00 POO rp 

+ J2 d P dqC m , p , q (t)\0) a \l p ,l q )\rrim 
m=0 Jo Jo 

where we have introduced the single- and two-photon dis- 
placed number states for the moving mirror 

|m(l)> 6 = e^ b+ - fe )|m>6, 

|m(2)) 6 = eW bt - b V) b , (3) 

with /?o = go I being the single-photon displacement 
quantity. We note that these displaced number states are 
eigenstates of the Hamiltonian 

-ffopc = hjj c a)a + hujM^b — Hgoa'aQr + b) (4) 

of the optomechanical cavity, as defined by the eigen- 
equation 

H p C \l) a \rh{l)) b = h(luj c + mcj M - l 2 v)\l) a \m(l)) b , (5) 

where v = g\juiM is the single-photon-state frequency 
shift. The energy-level structure limited in the zero-, 
one-, and two-photon subspaces of the optomechanical 
cavity is shown in Fig. [ljb). 

Based on Eqs. (JlJ, ©, and the Schrodinger equation, 
we obtain the following equations of motion for probabil- 
ity amplitudes 
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A m (t) = -i(mu] M - ±v)A m {t) - iV2^ / (m(2)\ b \h(l)) b B 7hk (t)dk, (6a) 

n=0 J 

oo oo /.oo 

An,feW = -«(A fc -^ + mw M )B m ,fe(i)-i\/2e^(" :i ( 1 )l f 'l"( 2 )) f ' A "W~^Zl / ("*( 1 )l6l n >6 c, n,p,Ai(t)tip, (6b) 

n=0 n=(K° 
oo 

C m>p , q (t) = -i(A p + A q + 



71=0 



Here we point out that the transition rates associated 
with photon scattering processes are determined by the 
Franck- Condon factors (m(2)| b |n(l)) b , (m(l)\ b \n(2)) b , 
(m(l)| b |n) b , and (m\ b \n( l)) b , which can be calculated 
based on the relation [23 



{m\ b e^^\n) b 



/ ^ e -^(-/3o)"- m ^ m (/3 2 ), n > m, 
" l (/3o), m > 



om — n j m—n i 
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where L s r {x) is the associated Laguerre polynomial. 



III. TWO-PHOTON SCATTERING SOLUTION 

Initially the cavity field is in the vacuum state and two 
photons are injected into the cavity. These two incident 
photons are prepared in a wave packet form outside the 
cavity. To facilitate analytic treatment by the Laplace 
transform, we examine the two-photon wave packet with 
Lorentzian spectrum. Without loss of generality, the ini- 
tial state of the mirror is assumed to be number state 
|no)t>. Once the solution in this case is found, the so- 
lution for general initial mirror states can be obtained 
accordingly by superposition. Explicitly, the initial con- 
dition is specified by A m (Q) = 0, B m k (0) = 0, and 



Cm,p,q(0) 



where Sj = i 2 
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(A p - 5 1 +ie)(A q - 5 2 + ie) 



+ S 1 «-> 6 2 



(8) 



ij — lu c (ojj is the resonant frequency) 
and e define the center detuning and spectral width of 
the two photons. The normalization constant J\f is 
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By using the Laplace transform, analytical solution of 
Eq. © subjected to the initial condition can be found 
(see Appendix,). In the long-time limit when the scat- 
tering is completed and the two photons exit the cavity, 
the solution is given by: A n0tm (oo) = 0, S„ , m ,fc(oo) = 0, 



and 

C„ , m , p>8 (oo) = Af[(Ci + C11 + Cm + Civ) 

+(A P <-> A q )]e-^ A " +A " +mUM)t (10) 

Here, we add the subscript no in A n0tm (oo), B nQtnii k(oo), 
and C„ 0imiPj g(oo) to mark the mirror's initial state \no} b . 
The four transition amplitude components C\, C\\, Cm, 
and Civ are given by 
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where we have introduced i*h-rv to denote products of 
Franck- Condon factors 

Fa = (m\ b \n(l)) b (h(l)\ b \n ) b , 

Fin = H 6 |fi(l))6<n(l)|6|n / >6<n / |6|r(l))6<r(l)|6|no)6, 

Fiv = (m| 6 |n(l)) b (n(l)| b |n'(2)) b (n'(2)| b |/(l)) 6 (Z(l)| 6 |n ) b . 

(12) 

In addition, the denominators in Eq. (| 1 1 1) are defined by 
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From Eq. (JT3J) we are able to determine the resonance 
conditions involved in the photon scattering process. 

We point out that the four amplitudes Ci_rv corre- 
spond to four different physical processes in the two- 
photon scattering (see Fig. [5]). The state transitions as- 
sociated with these processes can be identified by _Fn-iv- 
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FIG. 2: (Color online) (a-d) Four types of transitions from states |0)a|wo)ft to |0)a|'n}& of the optomechanical cavity in the 
two-photon scattering process, (e) The eigen-energies of these states involved in the four transition processes. 



For the term Cj, it describes a direct two-photon reflec- 
tion by the fixed end mirror without entering the cavity. 
Therefore the mirror state does not change, and the two 
reflected photons remain in the Lorentzian wave packet 
state [Fig. da)]. 

The term Cu corresponds to the one-photon scat- 
tering and one-photon reflection process, i.e., just one 
photon enters the cavity and the other photon is re- 
flected by the fixed end mirror. The entered photon in- 
duces the transition process \0) a \no)b — > |l) a |n(l))(, — > 
|0) a |m)h [Fig. H[b)], and such a process has the reso- 
nance conditions: n ujM + u>i — u c — v + uojm and 
u! p = (us c — v + tilum) — mu!M [Fig. E[e)], according to 
the poles of C u , i.e., Re(Mi) = and Re(M 1 - M 2 ) = 0. 

The third term Cm can be interpreted as a sequen- 
tial two-photon scattering process, in which the second 
photon enters the cavity after the first photon emitted 
out of the cavity. Such an interpretation is obtained 
from .Fin that the maximum cavity photon number is 
one. In this case, the system experiences the follow- 
ing transitions |0) a |ra o )& -> |l) a |Z(l)) b -> |0) a |ra')& -> 
|l)«|n(l)) 6 -> |0)«|m) 6 [Fig. [He)]. As depicted in 
Fig- Et e )i the photon excitation processes are governed 
by n uj M + u 2 = uj c - v + Iuj m [i.e., Re(M 4 - Ms) = 0] 
and ji'ujm = uj c — v + nuiM [Re (Mi — Mg) = 0], and 
the frequencies of the two resonant emitted photons are 
bj q = (oj c — v + IljJm) — n'ojM [Re(M4 — Ms) = 0] and 
u! p = (lu c — v + tiujm) — mu!M [Re (Mi) = 0]. 

The fourth term Civ corresponds to a 'genuine' two- 



photon process involving the interaction with two pho- 
tons inside the cavity. This is revealed in Fiy hav- 
ing the states with two cavity photons. The is a kind 
of two-photon cascade scattering in which the system 
experiences the transitions |0) a |no)b — > |l) |i(l))& — > 
|2)„|n'(2)) 6 -> |l),|n(l))i -> |0)«|m) 6 [Fig.^d)]. When 
the two photons are in the cavity, the mirror will expe- 
rience an energy shift —Ai>. This extra energy shift (and 
the phonon sidebands, as we will show below) provides 
the physical mechanism to create the photon blockade 
in the cavity. Due to the two-cavity-photon state is in- 
volved, the resonance conditions for the photon excita- 
tion processes are hquim + ^2 = u c — v + Iu>m [Re(M4 — 
M3) = 0] and (us c — v + Iujm) + = 2cj c — Av + n'uiM 
[Re (Ms — Mi) = 0]. In addition, the frequencies of the 
two emitted photons are uj q — (2u> c — Ai/+ti'ujm)~ (<^c — ^+ 
iicom) [Re(M6— Mi) = 0] andw p = (w c — v+nwfcf)— m% 
[Re (Mi) = 0]. These resonance conditions can also be 
seen from Fig. (He). 



The above analysis expatiated the physical picture for 
creation of two-photon correlation. In fact, after the scat- 
tering, the two photons will also entangle with the mirror. 
We can see this point by examining the scattering state of 
the system. We know from Eqs. and (|10p that, corre- 
sponding to the mirror's initial state |no)b; the long-time 
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state of the system is |0) a ® |$„ (oo)) with 

|$„ (oo)) = V / dp / dgC 7lD , m) p >g (oo)|m)h|l p ,l 9 ). 

(14 

Physically, the cavity is in a vacuum in the long-time 
limit, and hence it decouples with the environment and 
the mirror. The state |$„ (oo)) is an entangled state 
involving these outside fields and the mirror. We can un- 
derstand this entanglement from the dependence of the 
emitted photon modes on the mirror's final state \m)b- 
In the scattering process, the modes of the two emit- 
ted photons depend on the final state of the mirror [see 
Pig. He)]. 



IV. TWO-PHOTON JOINT SPECTRUM 

We know from the above section that the radiation- 
pressure coupling can induce photon correlation. We 
now show how to indicate qualitatively this correlation 
from the joint spectrum of the two scattered photons. 
Corresponding to the cases where the mirror is initially 



in a pure state \tp)i 



E 



:0 °™0 



n )b 



and a mixed 



state p b = J2n„=oP"o\ n o)b(no\b, the joint spectrum func- 
tions [6j are defined by 



S(A p ,A q )=J2 



m=0 



y ] c n Cn ,m,p,q{t 
n =0 



(15a) 



S(A p ,A q ) =J2Y, Pna\Cn ,m ,P,g(°°)l • ( 15b ) 
m— no— 

In Fig. [3l we plot the two-photon joint spectrum 
S(A P , A q ) as a function of the frequencies A p and A q 
when the initial state of the mirror is 1 0) 5 . We firstly con- 
sider the case of (3q <C 1 so that it is reasonable to approx- 
imately expand the probability amplitude C„ ,„ l;Pj9 (oo) 
up to the zero- and first-orders of /3q, i.e., using the ap- 
proximations 



(m\ b e^^\n) b 

(m| & e /3o(bt ~ b) |n)fc ~ $m,n + @oWn + 16 m>n+ i - VnS m , n -i) 

(16) 

in the zero- and first-order expansions, respectively. We 
note that the energy shift terms v and Av in Eq. (fT3|) 
are kept because v = 50 A) could be larger than "f c even 
for small 0q. Actually, the Kerr nonlinear energy shift 
!/(a'a) 2 is responsible for the generation of photon cor- 
relation. For the zero-order expansion, the effect of the 
sidebands is completely canceled, and the present op- 
tomechanical system reduces to a Kerr nonlinear cavity 
with the Kerr parameter v. When v > j c , two free pho- 
tons scattered by the Kerr nonlinear cavity will be cor- 
related @ . This can be seen from the joint spectrum of 
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FIG. 3: (Color online) Plot of the two-photon joint spectrum 
S(A P , A q ) for the ground state |0)& case. For a small /3o <C 1, 
we plot the spectrum by expanding the C no ,m, P ,q(oo) up to 
(a) the zero-order (u = 0.2) and (b) the first-order (/3o = 
0.4) of Pq. (c) For a large /3o, the joint spectrum includes 
many sidebands (/3o = 0.6). (d) The zoomed view of the 
peak with the center position A p /cjm — A,/wm = —0.36 
from the subfigure (c). Other parameters are: Jc/ujm = 0.1, 
e/ujM = 0.01, and 5i = 82 = —v- 



the two photons [Fig. [3]Ja)]. The two scattered photons 
are frequency anticorrelated with a probability concen- 
trated along the line parallel to A p + A q = 0. For larger 
fto, the phonon sidebands will be involved in the spec- 
trum. In the first-order expansion case, we can see two 
sideband peaks in the two directions of the joint spec- 
trum. With the increasing of /3q, more and more side- 
bands can be observed in the spectrum. In Fig.[3tc), we 
plot the joint spectrum in the single-photon strong cou- 
pling regime (3q = 0.6. From the spectrum, we can see 
that the pattern of these sideband peaks is concentrated 
along the line parallel to A p + A q = 0. This point can 
be seen clearly from Fig. Eld), which is a zoomed view 
of the peak with the center located at A p = A q = — v in 
Fig-EJc). We should emphasize that the joint spectrum is 
experimentally measurable by detecting the probability 
distribution of the two scattered photons. 



V. PHOTON BLOCKADE IN CAVITY 

In the previous section, we have studied the two- 
photon joint spectrum based on the long-time scatter- 
ing solution. Actually, we can quantitatively study the 
photon correlation by calculating the second-order corre- 
lation function for cavity photons. However, we now need 
to consider the transient dynamics rather than the long- 
time dynamics of the system, because there are no pho- 



6 




60 90 



FIG. 4: (Color online) Evolution of (a) the single-photon 
probability Pi (t) and (b) two-photon probability P2 (t) in the 
cavity for various go. (c) Equal-time second-order correlation 
function g^ (t) vs the scaled time um t. Other parameters are: 
■jc/iOM = 0.1, e/uiu = 0.01, c ng = 5 no ,o, and 5i = 62 = —v. 



tons in the cavity in the long-time limit. In the following, 
we calculate the transient dynamics of the system based 
on the solution given in the Appendix. Corresponding to 
these two cases of initial states \ip)h an d Pb, the proba- 
bilities for finding one and two photons in the cavity can 
be obtained as 



and 



Pi(t) 

W) 

Pl(t) 



E 



-n->0 



m— 

00 

E 



dk 



n =0 



n =0 



(17) 



J2Pn / dk \ R 



no— m— 



E P ™« E \ A no,m{t)\' 
no— m— 



(18) 



where A n0im (t) and B no ^ m _ k (t) are obtained by making 
the inverse Laplace transform on Eqs. (jA.101) and (|A.11I) . 



For the cavity field, the equal-time second-order corre- 
lation function is defined by [29[ 



2) /Vi - 



'(*) 



Tr[p(t)aia<aa 
(Tr[p(t)tfa\) 2 



(19) 



where is the density matrix of the total system at 
time t. In terms of Eqs. ©, ([17)1. and (fTgj) . Eq. (0 can 

be expressed as 



2P 2 {t) 



(20) 



Photon blockade effect corresponds to g^ 2 \t) <C 1 [30I ] . 
and the effect is measured by how small g^ 2 '(t) is. 

In what follows, we will study how the cavity pho- 
ton statistics depends on the radiation-pressure coupling 
strength g . In Fig.0] we plot the time evolution of Pi (t), 
P2H), and £f' 2 )(i) for various g at <5i = S- 2 = —u, which 
corresponds to a single photon resonance [see Fig. [T{b)]. 
Here we assume that the initial state of the mirror is 
\0)b, i.e., c„ = Sn ,o- We see that the probabilities Pi(t) 
and Pi{t) increase gradually from zero to a maximum 
value and then decrease to zero as time increases. At 
first glance, the maximum value of P2(t) is expected to 
decrease with the increase of go because the two-photon 
detuning —2v (i.e., the difference between the two-photon 
frequency shift —4v and two photon's total energy — 2v) 
is proportional to g 2 , and therefore the photon blockade 
effect becomes stronger as go increases. This trend is 
shown in the figure for go/^M from 0.01 to 0.6. How- 
ever, we see that the dependence of p2(t) on go is not 
monotonic due to the phonon sideband resonance effect. 
In fact, Fig. Ufc) shows that the photon blockade effect 
becomes weaker when go/^M = 1- We can explain this 
feature by the fact when go /(jJm = 1 , the two photons can 
induce the resonant transitions |0) a |0)& — !• |l) a |0(l))(, and 
|l>o|0(l)>6 -> |2) a |2(2)) 6j and so the maximum of P 2 (t) 
in this case is larger than that for gojiou = 0.6. 

The oscillating feature of the correlation function g^ 
as a function of go was found in Ref. [17| , which is based 
on an approximate steady state solution of a continu- 
ously driven system. In our case, our exact solution also 
indicates such a feature at transient times. This is shown 
in Fig. [S] in which the dependence of g^ (t) on go at a 
given time t p is plotted. Here t p is chosen when there are 
appreciable amount of photons inside the cavity. We can 
see clearly resonance peaks at specific values of go- To 
explain the resonance, we note that when the frequen- 
cies of two photons are 5\ = S 2 = — v used in the figure, 
resonant transitions |0) a |0)h — > |2) a |n(2))f, are allowed 
if — Af + uujm = —2v. This is the resonance condition 
allowing two photons to exist in the cavity at the same 
time. Specifically, the values of go associated with these 
resonance peaks are located at 



9o n 

= \/o> n = 0,l,2,- 



(21) 



which are consistent with the locations of the resonance 
peaks in Fig. [S] We know from Eq. (|2~T|) that, with the 
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FIG. 5: (Color online) Plot of the equal-time second-order 
correlation function g( 2 '(t p ) at time u)M~tp = 50 vs the scaled 
coupling strength go /cjm in the single-photon resonance case 
Si — &2 = —v. Other parameters are the same as those in 
Fig. H 



increasing of go, the sideband modulation peaks become 
more and more dense. 

The above discussions have considered the initial 
ground state of the mirror. It is interesting to ask 
how the g^ 2 \t) behaves when the mirror is in excited 
states initially. Suppose the initial state of the mir- 
ror is |n)b, then the system can undergo the transi- 
tions \0) a \n) b -> |l),|n(l)) 6 -> |2)„|m(2)) 6 . Therefore, 
the Franck- Condon factors (n\b\n(l)}b and (n(l)|j,|m(2))& 
are important to determine the magnitude of g( 2 \t). 
In Fig. |51 we illustrate this feature by considering a 
Fock state as an initial state (blue dot dashed 

curve). For the parameters used in this figure, the 
relatively small Franck-Condon factor (l(l)|i,|l(2));, < 
(0(l)|b|0(2))b leads to a suppression of i^i), and hence 
g^ 2 \t) can be substantially lower than that of the ini- 
tial ground state (red dashed curve). In Fig. [51 we have 
also plotted the g^ 2 \t) (black solid curve) for the initial 
thermal state pf{n = 1) = E^=o 2 ~ ( " 0+1) K)b("o|b (n 
being the average thermal phonon number). We see that 
the g( 2 \t) in this case is between those in the cases of 
1 0)b and 1 1) f> due to statistical mixture. 



FIG. 6: (Color online) The second-order correlation function 
g (t) vs the scaled time uiMt when the mirror's initial state 
is Fock state \l)b (blue dot dashed curve) and thermal state 
pl h (n — 1) (black solid curve). The ground state |0)b case (red 
dashed curve) is presented as a referrence. Other parameters 
are: <?o/Vm = 0.3, y e /u)M = 0.1, e/ujM = 0.01, and 5i — S2 = 
—v. 



two-photon correlation. In addition, from the poles of 
these scattering amplitudes, we have obtained the reso- 
nance conditions of these transitions. In particular, the 
final state of the two photons is generally correlated in 
frequency space because of the frequency sum appear- 
ing in the denominator of the amplitudes Cm and Civ- 
Based on the long-time solution, we have calculated the 
two-photon joint spectrum, which shows clear evidence 
for two-photon frequency anti-correlation. On the other 
hand, the transient dynamics of the equal-time second- 
order correlation function of the cavity photons has been 
calculated in order to address the photon blockade ef- 
fect. We have found that the photon blockade effect can 
be induced by the optomechanical coupling. Besides, the 
correlation function has also been found to exhibit reso- 
nance peaks as a function of the optomechanical coupling 
strength, and we have determined the peak positions in 
Eq. (|!?Tj) by a two-photon resonance condition. 

Note added. We notice a related paper that appeared 
recently [31 1. 



VI. CONCLUSION 

In conclusion, we have studied analytically the two- 
photon scattering in a cavity optomechanical system. 
Under the Wigner-Weisskopf framework, we have ob- 
tained the exact transient solution of the system with 
the Laplace transform method. On one hand, the long- 
time solution reveals the detailed physical transitions for 
the two photons in the scattering process. This could 
help us to understand the physical mechanism to induce 
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Appendix: Derivation of Eq. (|10p 

In this Appendix, we give a detailed derivation of the two-photon scattering solution in Eq. (|10p . By the Laplace 
transform f(s) = J °° f(t)e~ st dt, Eq. ^ becomes 



[8 + i(mu M -4v)]A m (s) = A m (p)-iV2£J2 / (M 2 )\b\n(l))bB n , k (s)dk, (A.la) 

oo oo ,.00 

[s + i(A k - v + mLo M )]B m ,k{ s ) = B m , k (0) ~ iV2£^(m(l)\ b \h(2)) b A n (s) - / (m(l)|b|n) 6 C n ,p i fe(s)(ip ! 

n=0 n=0^° 

(A.lb) 

oo 

[s + i(Ap + A q + mu M )}C m ,p,q{s) = C m ,p,q(0) - (m] &ln(l))t [B n ,p(s) + B n ,q(s)}, (A.lc) 



n=0 



where A m (0), B myk (0), and C mjPjg (0) are the initial conditions, which are given in Eq. flSJ. From Eqs. (|A.la|) . (|A.lc[) . 
and ((HI, we have 

£»(«) = r ^ .r JaT E / (m(2)| 6 |n(l)) 6J B n , fc ,( S )dfc', (A.2a) 



[s + i{A p + A q + mu M )\ V / 
Substitution of Eqs. (|A.2al) and (jA.2b|) into Eq. (|A.lb|) leads to 

s + — + i(A k -v + mcj M )J B m , k (s) 

y / 2^(m(l)| 6 |n(2)) 6 (n(2)| 6 |[(l)) 6 g 2 (m(l)|„|rc) 6 H,,|f(l)) 6 \ - 
^ I r- + i(nw M - 4i/)] + [a + i(A p + A k + mJ M )} J l * K)P 



Ln=0 \ L 



+2i^Af(m(l)\ b \n ) b ( — - ) — — -r +S 1 ^S 2 ), ( A.S) 

V[A fc + di + n WM - Hs + e)J(A fc - S 2 + it) 

where j c = 2-7r£ 2 is the cavity-field decay rate. We introduce a new variable F myk (s) by 

6 / \ 27r^A/"(m(l)|b|n )b / 1 , . . , , . „ 

[s + -f + i(A k -v + moj M )\ V [A k +h + n Ld M ~ H s + e )\{A k ~ ^2 + «e) 

Then the equation for F m ^ k (s) is obtained as 

(m(l)|b|n )b 



+ 5i o <5 2 F m , fc (s) 



[Afe + 5i + tiqWm - i{s + e)](Afe - <5 2 + ie) 

°° — / 2e 2 (m(l)| fc |n(2))b(n(2)|b|f(l))b | e(rh(l)\ b \n) b (n\ b \l(l)) b \ <f(l)| & | rao > & 

+ i{nu)M — 4f)] [s + z(A p + A fe + nwjw)] y [s + f + i(A p - c + Ium)] 



E 

;,n=o 



f 



+ <5i <-> 82) F lyP (s)dp + Si, (A.5) 



, [A p + (5i + lio^M — i(s + e)](A p — 8 2 + ie) 
where 

2 7c (m(l)|b|n(2))b(n(2)|b|[(l))b(/(l)|b|«o)fc 



[Si + 82 + n 0J M ~i{s + 2e)] /-^ 

Ln—O 



s + i(nu>M — 4z/)][#i - v + Iuj m — i(s + e + 
i7c<TO(l)|6|n)6(n|bK(l))6(Kl)|6|no)b 



[<5i + A fc + nwjvf — i(e + s)][5i — ^ + /cj m - «(s + e + 



(A.6) 
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By introducing a new variable x m (s) 



-l 



F m , k (s) = (jj- , f ( ^' fc| " o)fc -— T +S 1 ^5 2 ) [ft + x m (s)] , (A.7) 

we obtain the equation for x m (s) as follows: 

Xm{s)+ y lc(Mmm) b (m\ b m) b xi{c) = y 2 7e 2 (m(l)| b |n(2)) b (n(2)| b |[(l)) b (f(l)| b |n ) b 

+ (Si <-> S 2 ) , (A.8) 



[di — v + Iujm — i(s + e + ^ )] 



The solution of Eq. (|A.8I) can be found as 

2 7c 2 (m(l)| 6 |n(2)) 6 (fi(2)| 6 |r(l)) 6 <i(l)| 6 |no)6 



X rn I s 



) - E 



_ Q [s + 7c + i(«WM — 4i/)] [s + 2e + + £ 2 + n w M )][s + i(ru>JM — 4i/)] 

^- ; rT+5lO<S2). (A.9) 

[{s + e+^) + i{8 x -v + luj M )] J 
Then by Eqs. (|A~2a| . ([A~2bl) . (|A~4| . (TO) , and fO]) . we can obtain 

A , s) = 2V27rW 7c y> / (m(2)| b |n(l)) b (n(l)| b |n ) b | ^ ; > \ 

[s + 2e + i(Si + 62 +iioU)m)][s + 7c + i{muiM — 4z/)] ^ \[s + e+ ^ — v + uujm)] J 



5 f x = 2^(m(l)| b |n ) b / 1 1 + A ^ / 

m,M J [s + f + i(A k -v + mcJ M )} \ (Afc - <5 2 + ie) [a + e + i(A fc + ^ + now*)] 1 2 , 

_ g 47rzeAA 7c (m(l)| b |n(2)) b (n(2)| b |f(l)) b ([(l)| b |n ) b 



s + ^ + i(A k - v + mijj M )][s + 2e + i(Si + S 2 + ?i ^m)] 

+ (5l O 



[s + 7c + i{m)M - 4i/)] \ [s + e + If + i(<5i - 1/ + Icjm )] 
2 7 r^AA 7c (m(l)| b |n) b (n| b |/'(l)) b ([(l)| b |no) b 



E r„ 4. ic 



s + ^ + i(Afc - f + mw&r)][s + 2e + z(<5i + <5 2 + nowjf)] 

1 



[s + e + i(A k + <5i + nuj M )} [s + e + ^ + i{5i - v + lco M )} 



+ Si^5 2 ), (A.ll) 



Cm,p,g (^) 



A' 



[s + ?'(A P + A«j + mwAf )] 
1 



+ E 



X 



i 7e (TO| b |n(l)) b (n(l)| b |n ) b 
(Ap - Si + it) (A, - 5 2 + ic) "" v " ' ^ [a + f + i(A p - 1/ + nujif)] 



+ (5 X o <5 2 



[s + e + i(Ap + <5i + no^Af)] (A p - 5 2 + ie) 

7 2 (m| b |n(l)) b (n(l)| b |n') b (n'| b |f(l)) b ([(l)| b |n ) b 



E 



n^n'l— 



+ 2e + z(<5i + <S 2 + hjwm)][s + + *(Ap - v + hum)] 



[s + e + 3jf + «(<5i — v + Iuj m )] [s + e + «(A P + <5i + nVw)] 



+ 5i o 



E 

n.n'l— 



2 7e 2 (m| b |n(l)) b (n(l)| b |n , (2)) b (n , (2)| b |[(l)) b ([(l)| b |n ) b 

[s + 2e + + 5 2 + i!oWm)][s + 7c + i(n'u} M - 4i/)][s + ^ + «(A p - v + twjm)\ 



[s + e + ^ + i(5i —v + Iu>m)] 



Si O <5 2 



A p o A,. 



(A.12) 
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The transient solution of these probability amplitudes A„ 0)TO (t), B n0t7n ^(t), and C no ,m,p,q{t) can be obtained by the 
inverse Laplace transform. Here we add the subscript n in the transient solution to mark the mirror's initial state 
\no)b- The corresponding long-time solution is given in Eq. (fTO]). 
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